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Abstract 

Stochastic bistable systems whose stationary distributions belong to the (/-exponential family 
are investigated using two approaches: (i) the Langevin model subjected to additive and quadratic 
multiplicative noise, and (ii) the superstatistical model. Previously, the bistable Langevin model 
has been analyzed under linear multiplicative noise, whereas this paper reports on quadratic mul- 
tiplicative noise, which is more physically meaningful. The stationary distribution of the Langevin 
model under quadratic multiplicative noise, which agrees with that derived by the maximum Tsallis 
entropy method, is found to be qualitatively different from its counterpart under linear multiplica- 
tive noise. We also show that the stationary distribution of the superstatistical model is the same 
as that of the Langevin model, whereas their transient properties, described in terms of mean first 
passage times (MFPTs), are qualitatively different. 
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1. INTRODUCTION 



The exponential family given by 



p{x-(5) ■.= eM-m-m{x)], (1) 

is of great interest in physics, because of its relevance to many physical phenomena (/3 is 
inverse temperature, 9{(i) is a normalization term and H{x) is the Hamiltonian) . The expo- 
nential family includes the Gibbs measure, the canonical ensemble and the Gaussian distri- 
bution as special cases. The exponential family is important, since it appears as the limiting 
distribution of the central limit theorem, as distributions which maximize the Boltzmann- 
Gibbs-Shannon (BGS) entropy and also as stationary distributions of stochastic processes. 

In recent years, many investigations [l-7| have been made of physical phenomena which 
belong to the g-exponential family (such as the nonextensive canonical ensemble and the 
g- Gaussian distribution) js, 



p,(x;/3):=expJ-e(/3)-/3i/(x)], (2) 

where q is an entropic index and exp^(a;) is the g-exponential function [see Eq. ([7])]. The q- 
exponential family is a one parameter generalization of the conventional exponential family, 
and reduces to the exponential family as g — )■ 1. The g-exponential family can account for 
systems where physical quantities such as energy and entropy are not proportional to system 
size (nonextensive). As the conventional exponential family, the g-exponential family also 
plays an important role in stochastic processes. For example, particles under a quadratic 
potential driven by additive and linear multiplicative noise satisfy a g-Gaussian as their 
stationary distributions |lO |. 

In stochastic processes, bistable potentials are highly important in many fields includ- 
ing physics, electronics, biology and chemistry. It has been applied to chemical reactions, 
optical bistability, electric circuits, gene expression mechanisms and all the rest. Bistable 
stochastic processes can model switching dynamics of two-states systems under fluctuant 
environments. For bistable systems, multiplicative noise plays an essential role, and phe- 

12| and noise-enhanced stability 



nomena such as resonant activation 
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14| emerge 



only in the presence of multiplicative noise. In these studies, linear multiplicative noise has 
been exclusively assumed, since it is straightforward to analytically obtain its stationary 
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distribution. In this paper, on the other hand, we study quadratic multiphcative noise for 
a quartic bistable system. Considering a physical meaning of multiplicative noise, we show 
that quadratic multiplicative noise for the quartic potential is a straightforward extension 
of linear multiplicative noise for the quadratic potential. The quadratic multiplicative noise 
inherits an important property, namely that the noise intensity vanishes at stable sites, in 
the same way as intensity of linear multiplicative noise also vanishes at a stable site of the 
quadratic potential. In this paper, we investigate the stationary and transient properties 
of bistable quartic systems subjected to quadratic multiplicative noise. We show that the 
resulting stationary distribution belongs to the g-exponential family [Eq. ([2])], which agrees 
with a maximizer of the Tsallis entropy under constraints. Since much attention has been 
paid to the g-exponential aspect of physical phenomena in recent years, it is important to 
investigate bistable systems which belong to the g-exponential family. We first study sta- 
tionary distributions of the quadratic multiplicative case, comparing them with those of the 
linear multiplicative case for the quartic bistable potential. From these analyses, we show 
that the effects of multiplicative noise are different in each system. 



One of the alternative approaches yielding the g-exponential is superstatistics 15l-ll7|. 



Superstatistics can model or describe quasi-equilibrium systems, where environments fluc- 
tuate spatially and/or temporally. In recent years, Ref. [l8| derived the superstatistical 
Brownian motion from the viewpoint of mesoscopic nonequilibrium thermodynamics jl^. 
The concept of superstatistics has been extended to the path-integral 20| , and it was shown 
that the path-integral superstatistics include a financial model with stochastic volatility j2l| . 
These studies show that superstatistics is highly important for modeling and understanding 
of many real world phenomena. Regarding a relation to the g-exponential distributions. 



Ref. 



22] showed that Brownian particles moving in fluctuant environments satisfy the g- 



Gaussian as their stationary distributions. Therefore, we apply the superstatistical concept 
to the quartic bistable potential, and show that the stationary distributions belong to the 
g-exponential family in the small noise limit. 

In order to investigate transient properties of bistable stochastic processes in the g- 
exponential family, we calculate the mean first passage time (MFPT). We obtain the approx- 
imate analytic expression of MFPTs for the two cases of the Langevin model with quadratic 
multiplicative noise and superstatistics. We find that the g-dependence of MFPTs in the two 
cases is completely different, although their stationary distributions are the same. We also 
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see a formal nonextensive generalization in two MFPTs: exp(x) is replaced by exp2_q(x) in 
the quadratic multiplicative case and by exp^(x) in the superstatistical case. 

This paper is organized as follows: In Sec. [21 we first derive stationary distributions 
by the Langevin equation with quadratic multiplicative noise model and superstatistical 
approach. We next proceed to the calculation of MFPTs in Sec. [31 In Sec. [H we discuss 
the validity of our quadratic multiplicative noise from the viewpoint of an open quantum 
system. Effects of correlation between additive and multiplicative noise are also discussed. 
We conclude this paper in Sec. [51 

2. STATIONARY DISTRIBUTION 

2.1. Maximum Tsallis Entropy Principle 

In nonextensive statistics, many important distributions can be derived from the maxi- 
mum Tsallis entropy principle. The Tsallis entropy is given by 



dxP{xy 

S, := . (3) 

where q is an entropic index and Eq. (l3|) reduces to the BGS entropy in the limit as g — )■ 1. 



Adopting the optimal Lagrange multiplier (OLM) maximum entropy method 23| with the 
constraints, 

^ 1, (4) 

J dX^g{x)U{x) = Ug, (5) 

we obtain the distribution given by 

P{x) = exp^ [- ln2_, Z-f3 {U{x) - f/,)] . (6) 

Here U{x) denotes a potential, ^g{x) is the escort distribution defined by ^q{x) := 
P{xY / J dxP{xY, Z is a normalization constant and expg(x) denotes the g-exponential 
function defined by 

exp^(x) := [1 + (1 - (7) 
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(a) Quadratic potential with 
linear multiplicative noise 



(b) Quartic potential with 
quadratic multiplicative noise 



FIG. 1: (Color online) (a) A quadratic potential with linear multiplicative noise and (b) a quartic 
potential with quadratic multiplicative noise. Solid and dashed lines represent the potentials and 
multiplicative terms, respectively. 

where := max(a;, 0). Its inverse function, the g-logarithm, is defined by 

x^'i — 1 

Ing(x) := — (for x > 0). (8) 

With the use of these generalized functions, important nonextensive distributions can be 
expressed in a similar form to that of conventional ones. 

2.2. The Langevin Model Subjected to Quadratic Multiplicative Noise 

We here consider the following Langevin equation: 

- = f{x)+g{xm+vit), (9) 

where f{x) := —U'{x) and g{x) is a multiplicative term, ^{t) and ri{t) are white Gaussian 
noise with correlations: 

mas)) = 2M5{t-s), (10) 
{^{t)r^{s)) = 2A5{t-s), (11) 

mm = mvis)) = o, (12) 

where A and M represent intensity of additive and multiplicative noise, respectively. 
The Fokker-PIanck equation (in Stratonovich's sense) of Eq. iQ is given by 

^P(x, t) = -^F{x)P{x, t) + A_G(x)P(x, t), (13) 



where 



If g{x) and /(x) satisfy 



or identically 



Fix) = fix) + Mgix)g'ix), (14) 
G{x) = A + Mg{xf. (15) 

f(x) = -Kg{x)g'{x), (16) 



U{x) -Uo = ^gixy, (17) 



the stationary distributions are the g-exponential distributions as will be shown shortly [Eq. 



(k > is a proportional constant and Uq is the minimum of the potential) [10|. We 
study a simple quartic bistable potential given by 

U{x)-Uo = \{x'-l)\ (18) 
As a consequence, we consider the following multiplicative term: 

gix) = -^{x' - 1). (19) 



We next consider a physical meaning of the quadratic multiplicative noise of Eq. (]T9l) . 
For the case of stochastic processes with a quadratic potential and linear multiplicative noise 
(Fig. [U (a)), the stationary distributions are g-Gaussian. In this case, the nonextensivity 
is derived from the linear multiplicative noise [g{x) = x]. This linear multiplicative noise 
is natural in some cases, because the intensity of the multiplicative noise vanishes at stable 
positions and the intensity of the noise is greater when a particle is at more unstable states. 
Extending this property to the quartic bistable potential [Eq. ( 1T8|) ]. it may be natural to 
allow the intensity of multiplicative noise at the two stable sites to vanish. Because the 
multiplicative term defined by Eq. (fT9|) inherits this property (Fig. [I](b)), it is considered 
that Eq. f|T9|) is a straightforward extension of the linear multiplicative noise for quadratic 
potentials. In Sec. 14.11 a quantum interpretation for the multiplicative noise is discussed. 

We next calculate the stationary distribution. By calculating the stationary solution of 
the Fokker-Planck equation of Eq. f|T3|) . we obtain the following distribution: 
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Pst{x) oc exp {-V{x)) , 

with an effective potential V{x): 



(20) 



V{x) :-- 



du^J^ + lnG{x 



du 

K + M 



G{u] 
(M - K)g{u)g'{u) 



A + Mg{uY 

{U{x) - Uo) 



+ ln {A + MgixY) 



(21) 
(22) 

, - (23) 

2M V « 

Equations ( 120|) and (123|) are the g-exponential family [Eq. (|2])], since they can be rewritten 
as 



In M + 



PJx) 



exp. 



ln2-fZ-/3([/(x)-f/o) 



(24) 



where 



(25) 



k + M 

Z being a normalization constant. 

We plotted Eq. ( l20l) with Eq. ( I23l) in Fig. |2]with k = 2 {Z is evaluated using numerical 
integration). For comparison, we also show in Fig. |3] the stationary distributions for the 
linear multiplicative case [g{x) = x] with the quadratic potential given by Eq. ( IT8|) . for 
which an effective potential V{x) is given by 



For the case of quadratic multiplicative noise (Fig. |2]), the stationary distributions always 
exhibit bi-modality. We can see that the effect of M of quadratic multiplicative noise 
is different from that of linear one as shown in Figs. [2] and [3l For the case of linear 
multiplicative noise, larger M makes the probability density around x = higher. With 
M > 1, the stationary distribution for the linear multiplicative noise case is uni-modal (Fig. 
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FIG. 2: (Color online) Stationary distributions for quadratic multiplicative noise with each param- 
eter setting (k = 2). 




-2-1 I 2 3 -3 -2 -1 1 2 3 -3 -2 -1 



FIG. 3: (Color online) Stationary distributions for linear multiplicative noise with each parameter 
setting. 

[3]). For M < 1, the stationary distribution has two modals at a; = ±-\/l — M. This indicates 
that stable positions for the effective potential are not at x = ±1. On the other hand, 
the intensity of M behaves differently for quadratic multiplicative noise. In Fig. [2] (a), the 
density at stable sites (—1 and 1) is smaller for large M. Conversely, the density of stable 
sites (1 and —1) is higher for larger M in Fig. [2](c). From this result, we can see that effects 
of quadratic multiplicative noise depend on the intensity of the additive noise. 



2.3. Superstatistical Model 

We next describe bistable systems of the g-exponential family from the viewpoint of 



superstatistics. Superstatistics |l5l-ll7l] has been developed to describe quasi-equilibrium 
systems, where parameters describing environments fluctuate spatially and/or temporally, 
because many physical systems are inhomogeneous. Let h{P] Pq) be a PDF of /3, where 
/3q is a hyper-parameter of the distribution. In general cases, superstatistics is given by 
calculating the expectation of the Gibbs measure p{ei\(3) in terms of h{l3; f3o): 



8 



df3 p(£,|/3)M/3;/3o). 



(26) 



Equation ( !26l) defines the so-called "Type-B" superstatistics. By taking hidj-Po) as the 
prior distribution, superstatistics can be considered as a Bayesian framework [2J] {P{ei] f3o) 
corresponds to a posterior). For h{f3; f3o) = 6{f3 — f3o), Eq. ( 126|) reduces to the Gibbs measure. 

Superstatistics can also be applied to the description of Brownian particles, where the 
particles are in inhomogeneous environments. Ref. 22] has reported that superstatistical 
Brownian particles follow g-Gaussian distributions. 

Superstatistics assumes Brownian particles in inhomogeneous environments. The super- 
statistical Brownian model assumes the following Langevin equation, driven by additive 
Gaussian white noise: 



dx 
'dt 



(27) 
(28) 



where /(x) := — f/'(x) and a is noise intensity. We first calculate the stationary distribution 
in a given sub-system driven by additive noise, which is represented by 



Z{a) 



exp 



U{x)-U, 



a 



(29) 



with 



Zia) 



dx exp 



U{x) - Uo 



a 



■ia 



■ia 



where U{x) is given by Eq. ( IT8l) and Ia{z) is the modified Bessel function of the first kind {Uq 
is included inside the exponential). The stationary distribution of superstatistical Brownian 
particles is obtained by calculating the expectation of Eq. ( 129|) in terms of the PDF of a. 
Since it is difficult to analytically calculate the superstatistical distribution (Type-B) for Eq. 
( 129|) . we approximate Z{a) using the steepest descent method, i.e. 
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-/3a 



FIG. 4: (Color online) The x^-distributions with specific parameters (/3o and n). The mean and 
variance are given by /3o and 2/3Q/n, respectively. Distributions for larger n approach the 5- function. 



Z(a) 




(30) 



Here, xi and X2 are stable sites in double well potentials (xi < X2). In this case, two 
peaks are approximated using two Gaussian distributions (this is different from the steepest 
descent method used in MFPTs). Note that Eq. ( 130|) yields a reliable approximation for 
the case of low noise intensity {a is sufficiently small). 

In superstatistics, we often consider the fluctuation of /3a '■= 1/a instead of a itself. 
Superstatistics assumes that Pa is not a constant but fluctuates spatially and/or temporally. 
Here, we consider a temporal fluctuation. The spatial case can be calculated in a similar 
way, but an extra term is required for calculating the expectation |22i]. If the fluctuation 
is temporally macroscopic (fluctuates over a long time range), the corresponding stationary 
distribution is given by averaging over the distribution of (3a'- 



00 



Pst{x) (X I d(3ah{(3a;(3o) ^^^^^ ^ • (31) 







In superstatistics, three types of distribution are often used for h{f3a', f3o)- the x^- 
distribution (x^-superstatistics), the log- normal distribution (log-normal superstatistics) 
and the inverse x^-distribution (inverse x^-superstatistics). The support of these distri- 
butions is (0,00). The x^-distribution is particularly important, since x^ superstatistics 
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reproduces the g-exponential family. As a consequence, we assume that /3a is sampled from 



the x^-distribution with n degrees of freedom (Fig. H]) 15|, |l6|. The assumption of the 
X^-distribution can be understood as follows: suppose that there are n independent random 
processes Yi, which are sampled from Gaussian distributions, behind If Pa is realized 
as a sum of squared Yi, then Pa is a random variable of a ^^-distribution with n degrees of 
freedom. The PDF of Pa is given by 



h{Pa;Po) 



\ "/2 



exp 



nPc 



(32) 



T{n/2) \2P0J "° V 2/3o, 

where /3o is a hyper-parameter, the average of Pa- The variance of the x^-distribution is 
given by 2Pl/n. By substituting Eq. ( 132|) in Eq. ( 13T1) . we obtain 



Pjx) = 



1 + ^ {U{x) - Uo) 
n 



n + l 
2 



(33) 



where is a normalization constant. Equation ( 133|) is a reliable solution when /3o and n 
are sufficiently large. Equation fl55]) can be represented by 



Psi(a;) = exp~ 



QB 



Pl 



Pb {U{x) - Uo) 
n + 3 



(34) 
(35) 



n " ~ n+l 
Equation ( l34|) is equivalent to Eq. ( l24|) . which indicates that superstatistical bistable 
stochastic processes also follow stationary distributions belonging to the g-exponential fam- 

iiy. 

The above calculation is carried out using Type-B superstatistics. We may alternatively 
calculate stationary distributions using Type-A superstatistics, where the factor of Z{l/Pa) 
in Eq. f l3T|) is neglected. This case also yields the g-exponential family given by Eq. 
but with 

PA = PoZf-\ U = '^. (36) 
where Za is a normalizing constant. 



n 
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(a) (b) 




M q 



FIG. 5: (Color online) MFPTs of the Langevin model with quadratic multiplicative noise {n = 2) 
as functions of (a) M and (b) q. 
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FIG. 6: (Color online) The effective potentials V{x) — Vq with k = 2. 
3. MEAN FIRST PASSAGE TIME 



3.1. The Langevin Model Subjected to Quadratic Multiplicative Noise 

In order to study transient properties of bistable systems, we calculate the MFPT, which 
is an average of the first passage time (FPT). The FPT is time required for a particle to 
arrive at one stable site from the other stable site. 



According to Fox's Ansatz 25|, an analytic expression for the MFPT of r from Xi to X2 
is expressed by 



X2 



where the Fokker-Plank equation of Pst{x) is given by Eq. ( !T3|) . Let r be a MFPT for the 
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quadratic multiplicative noise case. A Kramers-like formula 26| is obtained by applying the 
steepest descent method to Eq. fl37|) : 



2-K 



G{Xm)VV"{Xi) \ V"{Xm)\ 

where unstable site {xm = 0) and Vq := V{xi) [V{x) is given in Eq. (123|) ]. By 

substituting in Eq. fj38|) . we obtain the following expression: 



27r 



G{x„,)^V"{x^)\V"{xJ\ 

27f 



2M 

1 + — {U{Xm) - Uo 



K + M 

2M 



kA 



exp2_ 



/3o (f/(a:„) - f/o) 



(39) 



G'(a;^)VV^"(xi)|l^"(x^)| 
where (3q and g are defined in Eq. (12^ . and f/o = [/(xi). 

Equation (!39|) is a good approximation for V{xm) — Vq ^ 1. In the limit as M — )■ 0, Eq. 
fl39|l reduces to that for the conventional case. The last part of Eq. fl39|) is similar to that 
of the Kramers-time [Eq. (14T|) ]. with exp(x) replaced by exp2_q(x). We plotted Eq. (!39|) 
in Fig. [5] (k = 2). Fig [5] (a) and (b) show MFPTs as functions of M and g, respectively. 
In Fig. [5] (a), the MFPT decreases as the noise intensity increases. In Fig. [5] (b), it is also 
decreasing as a function of q. If n is constant, q depends only on the multiplicative noise 
intensity M. Because M is an increasing function for 1 < g < 3, Fig. [5] (b) decreases as q 
increases. 

For the case of linear multiplicative noise, a Kramers-like equation is obtained over a 
very narrow parameter range. As can be seen in Fig. [3l the effective potentials of stationary 
distributions cannot be well approximated with the steepest descent method when M is 
large. 



3.2. Superstatistical Model 

We next calculate the MFPT resulting from the superstatistical description. We assume 
that time scale of fluctuations of Pa is macroscopic, i.e. f3a does not change during each es- 
cape event. Under this assumption, the MFPT for the superstatistical case can be calculated 
by taking the expectation in terms of h{f3a', /3o): 

POO 

T= / d(3aTi(3a)h{(3a,(3o), (40) 
Jo 
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FIG. 7: (Color online) MFPTs of the superstatistical model as functions of (a) and (b) qj^. For 
comparison, we plotted the Kramers-time in (a). 



where T(/3a) is the MFPT of the Brownian model given by Eq. (1271) . 

27r 



exp [/3c, {U{x^) - f/o)] • 



(41) 



Note that Eq. f HTj) is the Kramers time and is valid for sufficiently large fi^. Since the 
Kramers time is valid for large /3q,, we obtain the following result for sufficiently large n and 
/So: 



27r 



(42) 



with n/{2[3Q) > U{xm) — Uq. We see that Eq. fj42|) is the same as Eq. pij) . except that exp 
is replaced by exp^. For the quadratic multiplicative case, exp is replaced by exp2_g. We 
see the dual relation g ^ 2 — g which often appears in nonextensive statistics. We plotted 
Eq. in Fig. [71 Figures [7] (a) and (b) show MFPTs as functions of := l//3o and 

Qa = {it- + 2)/r;,, respectively. When n — t- oo, r reduces to the Kramers-time. Fig [7] (b) 
shows that the MFPT is increasing as a function of qa (and also increasing as a function of 
Qb = (3$A — l)/(^4 + !))• As shown in Sec. 13. the MFPT for the quadratic multiplicative 
case is decreasing as a function of g. It is interesting to see that effects of the entropic index 
on the MFPTs for both cases are inversely related, although their stationary distribution 
dependence on the entropic index agrees. 

We have shown that the g-dependence of the MFPT in the Langevin model driven by 
quadratic multiplicative noise is different from that in the superstatistical model, which will 
be qualitatively explained as follows. In the Langevin model with quadratic multiplicative 
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noise, larger M with fixed A yields a smaller MFPT (Fig. [5]). Since larger M leads to larger 
q [Eq. f l25|) ]. the MFPT is decreasing as a function of q. On the contrary, the MFPT in 
the superstatistical model is expressed as superposition of MFPTs of given subsystems, in 
which smaller a [i.e. larger leads to larger r(/3Q,). The x^-distribution with smaller n 
(larger q^ [Eq. (1361) ]) has higher magnitude in large- /3a regions. Thus the MFPT in the 
superstatistical model is increasing as a function of g^i, which is opposite to that in the 
Langevin model. However, if we adopt an effective potential V{x) derived from the PDF of 
Eq. dMD 



V{x) 



n 



1 



In 



1 + ^ mx) 

n 



and substitute Eq. ( I43l) to Eq. ( 138|) . the exponential part of Eq. ( l38i) is given by: 

'n + 1 



exp [V{x.m) - Vo] = exp2^~ 



IB 



n 



-MUiXm) - Uo) 



(43) 



(44) 



The g-dependence of Eq. (jH]) is similar to that of the Langevin model with quadratic 
multiplicative noise mentioned above [Eq. ( l39l) ]. It is stressed that in order to study the 
MFPT (a typical dynamical quantity) in the superstatistics, it is necessary to take an average 
of T(/3a) over /3q, because the characteristic time of fluctuations in Pa is much slower than 
that of the MFPT. 



3.3. Wall Height Dependence 

We study the dependence of the MFPT on wall height {i.e. U (xm) — Uq). Let us consider 
the following functions: 

U{x) -Uo = a(x2 - 1)^ (45) 

g{x) = ^-{x' - 1), (46) 

where the coefficient a expresses the wall height U{xm) — Uq. We calculate the a-dependence 
of the MFPTs for the two statistics while keeping the noise terms unchanged. Eq. (fT7|l with 
Eq. (HE]) is identical to Eq. (jlSl) with k = 8a. 

Figures [H] (a) and (b) show the wall height dependence of the MFPT for the Langevin and 
superstatistical models, respectively. Since Fig. [H] (a) and (b) are log-plots, the exponential 
dependence on a is indicated by a straight line. Since the Kramers-time is dominated by an 
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(a) 



10' 



10= 



Kramers time^ 

<-'M = 0.1 



M= 1 




a 

FIG. 8: (Color online) MFPTs as a function of the wall height a. (a) the Langevin model with 
quadratic multiplicative noise and (b) the superstatistical model. In both figures, the intensity of 
additive noise is ^ = = 0.1. 

exponential of the wall height, a linear relationship is indicated in both figures (Although 
the Kramers-time depends on curvature of the potentials, its effect is smaller in comparison 
to the exponential part). Figure [S] (a) show that the log-scaled MFPT as a function of 
a exhibits a linear relation in the Langevin model as the Kramers time. A wall height 
dominant part of the MFPT [Eq. (139|1 ] for the quadratic multiplicative case is calculated as 
follows: 

A,(f/(xJ-t/o) 



exp2_ 



Sa + M 



(47) 



Eq. f H7|) shows the reason for the linear relation in the log-plots. 

We note that the a-dependence of the MFPT for the superstatistical case in Fig. [S] (b) 
is different from that of the Langevin model in Fig. |8] (a). For n = 30, the MFPT grows 
super-exponentially as a function of a. Equation f H2|) is dominated by the g-exponential 
function, and qa is large for small n {q^ = 1.07 for n = 30). Although the stationary aspects 
are the same in some limits, their MFPTs as functions of the wall height are qualitatively 
different. 
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4. DISCUSSION 



4.1. Open Quantum System 

A microscopic origin for multiplicative noise was discussed for an open quantum system 
It was assumed that a harmonic oscillator is coupled with harmonic oscillator bath by 
the interaction, 



where x and qj denote the coordinates of the system and bath, respectively, Cj is the coupling 
strength and f{x) the coupling function. It has been shown that the noise term for a quantum 
Langevin equation in the Markovian limit is given by ^] 



which yields additive [g{x) = 1] and multiplicative noise [g{x) = x] for f{x) = x and 
/(x) = x^/2, respectively. It is possible to extend the above discussion to an open bistable 
quantum system which is in a stable state of either x = x* = — lor+1. The coupling 
interaction may be expressed as 



for which we obtain a multiplicative noise of g{x) = x — x* for f{x — x*) = {x — x*y /2 in 

the relevant quantum Langevin equation. This implies that the multiplicative noise is the 

most pronounced at the stable state and that it is consistent with Eq. flT9|) because 

1 

g{x) = -{x^ — 1) ^ (x — X*) (for x — x* <^ 1). 
Thus our choice of g{x) given by Eq. ( IT9l) is consistent with the open quantum approach 



4.2. Effects of Correlation Between Additive and Multiplicative Noise 

In many researches on bistable systems using Langevin models, the correlation between 
additive and multiplicative noise is considered. In such cases, the correlation function be- 
tween additive and multiplicative noise is given by 



Hi = - ^^cjqjfix), 



j 



gix)v{t) = {f'ix)Ht), 




27|. 



mv{s)) = {vims)) = 2XVAM6{t - s), 



(48) 
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FIG. 9: (Color online) Stationary distributions for correlated additive and quadratic multiplicative 
noise with each parameter setting {k = 2). 




-3-2-1 12^-3-2-1 123 



FIG. 10: (Color online) Stationary distributions for correlated additive and linear multiplicative 
noise with each parameter setting. Results for A < are given by a symmetric change of those for 
A > with respect to the x = axis. 

instead of Eq. ( fT2l) . In Eq. ( I48l) . A represents the correlation intensity (|A| < 1). The 
Fokker-Planck equation [Eq. (fT3|) ] is given by 

F{x) = fix) + y/AMXg'{x) + Mg{x)g'{x), 

G{x) =A + 2\VAMg{x) + Mg{xf. 

For the linear and quadratic multiplicative cases, the stationary distributions can be obtained 
analytically with Eq. f l2T|) (the expressions are not shown here). We plotted the stationary 
distribution for correlated cases in Fig. |9] (quadratic multiplicative case) and Fig. [10] (linear 
multiplicative case) with k = 2. For the quadratic multiplicative case (Fig. [9]), stationary 
distributions are symmetric for every A. On the other hand, the linear multiplicative case 
results in asymmetric distributions, as shown. Furthermore, for a negative correlation region 
(A < 0), the stationary distributions for the linear multiplicative case are only line-symmetric 
distributions with respect to x = 0. In contrast, the quadratic multiplicative noise cases yield 
different shapes for negative A. Furthermore, it would be interesting to apply our method 
to stochastic systems with an asymmetric bistable potential which have been extensively 
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studied in recent years |28|, |29 |. 



5. CONCLUDING REMARKS 

In this paper, we have investigated properties of stochastic bistable systems described by 
the g-exponential family as given by Eq. ([2]), using two approaches: the Langevin model 
driven by quadratic multiplicative noise and superstatistics. We have pointed out that 
quadratic multiplicative noise is more physically appropriate for bistable quartic potential 
systems than the linear version and that it is consistent with the result of an open quantum 
system {2^. Properties of the stationary distribution for quadratic multiplicative noise are 
quite different from those for the linear one. We have shown that the MFPTs in the Langevin 
model and the superstatistical model are qualitatively different. 

Langevin equations (especially the bistable cases) are widely used in biological systems, 
where environments fluctuate temporally and / or spatially. Recently, superstatistics has been 
applied to a cancer survival modeling |30|, and it was shown that the superstatistical de- 
scription succeeded in modeling cancer survival in an accurate manner. This result indicates 
the importance of the superstatistical modeling for biological mechanisms. Since we only 
calculate here the superstatistical model for macroscopic fluctuations, it may be desirable 
to investigate the properties of more general cases. These are left as our future studies. 
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